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ABSTRACT 



Models of the Universe like the Concordance Model today used to interpret cosmological observations give expectation values for 
many cosmological observables so accurate that frequently peoples speak of Precision Cosmology. The quoted accuracies however 
do not include the effects of the priors used in optimizing the Model nor allow to evaluate the confidence one can attach to the Model. 
We suggest an estimator of the Confidence Level for Models and the accuracies of the expectation values of the Model observables. 
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1. Introduction 



After the discovery of the Cosmic Microwave Background 
(CMB) in 1964 (IPenzias and Wilson 19641 ). cosmological ob- 
servation went through a very rapid expansion. Measurements 
confirmed to very high degree of accuracy the expected first 
order (planckian spectrum, isotropic distribution, absence of 
polarization) and second order (angular power spectra of tiny 
temperature anisotropies and even smaller level of residual po- 
larization) features of the CMB (1COBE 20001 IWMAP 20101 . 
Complementary information then arrived from no-CMB chan- 
nels of information, namely: i)peculiarities in the rotation curves 
of galaxies and clusters of galaxies, (which suggest the existence 
of "invisible" or "dark" matter), ii)discovery of BAO (Barion 
Acoustic Oscillations) in the spatial distribution of the luminous 
matter, similar to the oscillations observed in the CMB angular 
power spectra, iii)Supernovae la used to get the distance of ob- 
jects at very large Z, iv)N-body simulations of the formation of 
matter condensations from isotropic distribution of matter parti- 
cles (1LSS 20101 1. 

The whole set of CMB and no-CMB information, besides 
confirming the general layout of the Standard Big Bang Model 
dPeebles 1 993). suggests that the Universe is probably suffering 
re-acceleration and its geometry substantially flat (euclidean). To 
close the Universe it was therefore assumed the existence, beside 
barionic and dark matter, of a third component of the matter- 
energy mixture, Dark Energy. Moreover besides the Hubble con- 
stant H , the ratio between the primordial abundances of he- 
lium and hydrogen and the Universe average matter density O 
(in units of the critical density p c = 3H}j{%nG)), new observ- 
ables were introduced: the abundances (still in unit of critical 
density) of barionic matter D./,, dark matter Q^ m and dark energy 
£2a, the optical thickness r of the Universe at reionization (when 
at 10 < Zi on < 1300 stable matter condensation emerged from 
the primordial uniform matter distribution), the amplitude and 
spectral index ( A s and n s ) of the fluctuations (seed for the birth 
of the above condensations), and erg an indicator of the fluctua- 
tions of the galactic matter distribution. 



To reconcile these new data and the old Big Bang model 
the Concordance Model or A-CDM Model gradually emerged 
(IQstriker and Steinhardt 19951 IKowalski et al. 20081 It has six 
independent parameters, representative of six observable chosen 
among the observables listed above, usually H„, Q;,, Qj m , r, A s 
and n s , from which other observable (e.g age of the Universe 
T„niv - l/Ho> critical density p c of matter-energy, Dark Energy 
density, Q.^ (in unit of p c ), reionization red shift Z io „ and <x 8 ) can 
be derived. 

The observables-parameters and the observational data are 
entangled: different sets of data are produced by different pa- 
rameter combinations. Disentangling the effects of the various 
parameters and getting the parameters values is extremely hard. 
A way out is obtained fitting the Model to the full set of CMB 
and no-CMB cosmological data and getting the parameter val- 
ues by best fit methods. This is done by Montecarlo meth- 
ods using Markov chains to implement the stochastic procedure 
which drives Montecarlo calculations dLandau and BindeF 2009) 
with additional a priori assumptions, e.g. the Universe geometry 
(Q = Q;, + D-dm + Qa = 1 if the Universe is fiat) and/or intervals 
of allowed variability for the parameters. The parameter values 
which give the best fit are interpreted as expectation values of 
the associated observables. 

This procedure, which forces the model toward pre- 
ferred solutions, is now well established. Frequently im- 
proved expectation values are published. They are ob- 
tained adding new observational data to the preexisting 
data base (e.glKomatsu et al. 20091 IBrown et al.~20"09l 
Komatsu et al. 20101) . Recently published values 
dKomatsu et al. 20T0b of free and derived observable - pa- 
rameters of the Concordance Model are shown in Table Q] and 
Table [2] The accuracy of the expectation values is so high 
that observers sometimes speak of Precision Cosmology (e.g. 
IPrimack 20051 IBridle et al."2003l . 

In principle all the observables are measurable. Agreement 
between measured and expected values is a proof that the model 
described by the free parameters is a good one. Unfortunately for 
five of the six free observables - parameters of the Concordance 
Model measurements are very poor or not yet available. And 
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no significant improvement is expected in the near future. For 
these observable in literature there are only large intervals in- 
side which the measured values are expected. The same inter- 
vals are usually assumed as the maximum variability range of 
the parameters used in Monte Carlo studies of the Concordance 
Model (B rown et al. 20091 . They are shown in TableQ] The only 
exception is the Hubble constant for which almost direct mea- 
surements are now very accurate (IRiess et al. 20091 1. 

2. Expected and measured values 

Measured values M' and their uncertainties 5 M , are ex- 
tracted using classical statistical methods (see for instance 
Bevington and Robinson 1992) ) from sets of repeated measure- 
ments of a parameter i or from combinations of measured quan- 
tities directly linked to i. The distribution of the measured values 

around their average M' is usually gaussian. 

Expected values E' v are the values of the model parameters 
which gives the best fit of the Model to the full set of observa- 
tions assuming a priori conditions or priors. The uncertainties 
6e v of E v are the widths of the distributions (not always gaus- 
sian) which encompass 68 % of the parameter values obtained 
repeating the calculation with a random initial choice inside a 
priori chosen intervals of free parameter variabilities. 

In doing it, sometimes unintentionally, a transition from clas- 
sical to bayesian statistics occurred (Stanford 2003). Bayesian 
probabilities, no longer precisely defined as frequencies or ra- 
tios of measured quantities, give the (0-1) confidence level we 
associate to an event occurrence and are obtained by a priori 
assumptions. The Bayes Theorem allows to evaluate the effects 
new observations have on preexisting knowledge: 

P E (H) = [P(H)/P(E)]P H (E) (1) 

where 

i) H is a set of measured values of a parameter and P(H) its 
probability distribution, derived for instance from previous ob- 
servations, assumed ab initio, i.e. before new observations are 
made; 

ii) E is a new set of measurements of the same parameter and 
P(E) its distribution assumed a priori, or prior, without preex- 
isting information; 

iii) Ph(E) is the conditional or "direct" probability of getting 
E given H; 

iv) Pe(H) is the a posteriori conditional probability of get- 
ting H given E, or "inverse" probability. 

The ratio < [P(H)/P(E)] < 1 is the Confidence level C/(£) we 
can associate to the prior used to improve our knowledge. 

For each parameter i of a multi parameter model we can write 

P'(H) = PL(XL< e ) (2) 

P\E) = PJKM (3) 

where P' me (X' me , cr l me ) and P' ex (X' ex , o-' ex ) are the probabilities dis- 
tributions of the measured (pedix me ) and expected (pedix ex ) 
values X', with dispersion <x'. Assuming gaussian distributions 
the Confidence level of the expected value of parameter i can be 
written: 



where AX' = \X< - X'\. 

For a model not too far from the real world 
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(6) 



-f exp{-[(AX' m Jcr' me ) z - (AX'Ja' ex ) z ]/2} 



The resulting Confidence Levels of the expected values of 
the free parameters of the Concordance Model are shown in the 
last column of Table Q] For the whole Model 

C, -0.4 

a value which leaves room for other models to be considered, 
obtained for instance relaxing priors, e.g. + Q f + Qa = 1. 



3. Discussion 

The numerical values of C\ and C\ presented above have been 
obtained for the A-CDM Concordance Model of the Universe. 
The difference between measured values and expectation values 
of the observable - parameters however hold for other Models. 
For each model we can evaluate expected value and Confidence 
Level C, of each observable and C/ for the whole Model. The un- 
certainties of the calculated Expectation values can be assumed 
as uncertainties on the observable real values only when the 
model priors are such that cr me - cr ex . 

If not, systematic uncertainties introduced by the priors have 
to be added to the uncertainties of the observable. 

A way to validate priors and model is therefore repeating di- 
rect (or almost direct) measurements of the parameters until their 
accuracies becomes comparable to the accuracy of the expected 
values. 

Until C; < 1 models based on different priors cannot be ex- 
cluded. 



4. Conclusion 

The use of Montecarlo methods and Bayesian Statistics to an- 
alyze the enormous quantity of data of cosmological interest 
which are continously poured by ground and space observa- 
tions is almost unavoidable. However Montecarlo and Bayesian 
Methods are based on a priori assumption whose statistical 
weight should be added, but rarely is added, to the the quoted 
accuracies of the parameter expected values. 

Peoples who currently use these methods are aware of that 
and warning have been put forward (e.g. IBridle et al. 2 003). 
Unfortunately general public and professionals not involved in 
cosmological observations may be unaware of it, misinterpret 
the results of model simulations and attribute weights above their 
real values to models. Forgetting it may stop or reduce support 
to studies of other models based on different priors still possible 
and not yet excluded by observation. 

Similar situations occur in other fields of pure and applied 
research (e.g. unification of fundamental forces, string theories, 
elementary particle models, models of climate evolution and so 
on). To avoid misunderstandings and preserve possibilities of 
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Table 1. ACDM-Concordance Model: Expectation values, Measured Values and Confidence Level of Model Parameters - (adapted 
from lKomatsu et al. 20101 IBrown et al.~2 009 ) 



Parameter < E > M C ; 

Expec. value Meas. value Conf. Level 

Hubble Constant (km/sec Mpc) H„ 70.4+m 74.2 ± 3.6 0.38 

Barionic Matter Density Q. b 0.0456 ± 0.0016 0.005-0.1 < 10~ 2 

Dark Matter density tl dm 0.227 ± 0.0014 0.006- 1 < 10~ 2 

Optical thickness at reionization r 0.087 + 0.0014 0.01 - 0.80 < 10~ 2 

Scalar fluctuations Amplitude A s (2.441!°;°** 10~ 9 ) ? ? 

Scalar Spectral index n s 0.963 ± 0.012 0.5-1.5 < 2 10~ 2 



Tabl e 2. ACDM -Concordance Model: Expectation values of Derived Parameters - (adapted from IKomatsu et al. 20101 
IBrown et al. 20091 ) 



Parameter Expected value 

Dark Energy Density Q A 0.728!° ^ 

Reionization Red Shift Z,,„, 10.4 ±1.2 

galactic fluctuations amplitude <xg 0.809 ± 0.024 

Universe Age (years) t (13.75 + 0.11) 10 9 



pursuing alternatives lines of research, this point of view must Riess, A.G. et al. (2009) ApJ699, 539 

be transmitted to the bodies who support research activities and http://plato.stanford.edu/entries/ bayes-theorem/| and references therein 

to the general public. http://wmap.gsfc. nasa.gov/| and references therein 
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